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Abstract 



The renormalised quark Dyson-Schwinger equation is studied in the limit of 



Q.C the renormalised current heavy quark mass iriFt ~> oo. We are particularly 

J~£ • interested in the analytic pole structure of the heavy quark propagator in the 

, complex momentum plane. Approximations in which the quark-gluon vertex 

£S) \ is modelled by either the bare vertex or the Ball-Chiu Ansatz, and the Landau 

gauge gluon propagator takes either a gaussian form or a gaussian form with 

►> ■ an ultraviolet asymptotic tail are used. 
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I. INTRODUCTION 

The solution of approximate Dyson-Schwinger equations (DSE) has proved to be an 
effective means for modelling quark propagators in hadronic physics |L) . Recent calculations 
within the genre of models which we shall refer to as the DSE technique include those 
of the light hadron spectrum 0,[3] and of electromagnetic form factors of the pion and 
kaon Q . Although no rigorous proof exists, it is the philosophy of the DSE technique that 
one possible signal of confinement in QCD should be the absence of timelike poles in the 
quark propagator ||. It has furthermore been conjectured that the propagator S(p) could 
be an entire function in the complex p 2 plane H0]. Such a scenario would, for instance, 
avoid certain unpleasant consequences which can result when modelling mesons via the 
Bethe-Salpeter equation which samples the quark propagator over a region of the complex 
plane. 

Determining the analytic structure of fermion propagators in QCD ||£| or other confining 
theories [|K],|11[] by the direct solution of model DSEs is not easy. It appears that the 



pole or branch cut structure obtained in any particular model is heavily dependent on the 
approximations employed. In general, two aspects of the quark DSE must be approximated: 
the quark-gluon vertex, and the gluon propagator. In this paper, we shall look at both 
these aspects within the heavy quark sector. In existing numerical studies in the light quark 
sector which produce propagators with conjugate singularities |8]||, the quark gluon vertex 
has usually been approximated by the bare vertex (the so called rainbow approximation). In 
ref. 0, however, it was shown that an entire function propagator can be obtained if the vertex 
function is modelled by a more sophisticated form respecting the Ward-Takahashi identity. 
This suggests that it is worthwhile exploring the importance of accurately modelling the 
quark-gluon vertex, as well as the gluon propagator, when studying the analytic structure 
of the quark propagator. 

In a recent development, the DSE technique has been extended to the realm of heavy 
quarks ]|12| , |I"3"1 in a way inspired by heavy quark effective theory (HQET) flUf . The purpose 



of this exercise was twofold. Firstly, if one acknowledges the success of the DSE technique 
in the light quark sector, it is clear that the dynamics of confined particles is driven by 
non-perturbative dynamical self dressing. In HQET, non-perturbative self dressing and the 
detailed analytic structure of the heavy quark propagator are largely ignored. It is important 
to know whether this is justified, or whether the successes of HQET are purely fortuitous. 
Secondly, one has the hope that an accurate determination of the heavy quark propagator 
will eventually prove useful for building phenomenological models of heavy quark hadrons. 
In ref. |12| a preliminary attempt is made to calculate the spectrum of heavy quark- 
light antiquark mesons by using the combination of rainbow DSE and ladder Bethe-Salpeter 
equation (BSE). It is found that, within the limitations of the model, the pole structure of the 
heavy quark propagator prevents solution of the meson BSE. This is clearly a shortcoming 



of the approximations involved. In ref. [13] the heavy quark DSE is examined from the 



point of view of the gauge technique. This is essentially an improvement on the rainbow 
approximation to the quark-gluon vertex which is designed to respect the Ward-Takahashi 
identity. An alternative approach, and one which we follow in this paper, is to replace the 



bare vertex Ansatz with the Ball-Chiu vertex Ansatz [15]. We shall see that in the heavy 



fermion limit the gauge technique and the Ball-Chiu vertex are equivalent. 



Regarding the gluon propagator, our treatment differs from ref. ||T3|] in that it is princi- 
pally numerical, enabling us to concentrate on a more realistic class of model gluon prop- 
agators. Specifically we study the simple gaussian Ansatz employed in ref. [T^], designed 



to model the infrared enhanced behaviour of the gluon propagator, and a model proposed 
by Frank and Roberts (FR) (16| which includes both an enhanced infrared behaviour and 
the known asymptotic ultraviolet behaviour. In order to deal with the FR propagator, it 
has been necessary to formulate a properly renormalised version of the heavy quark DSE of 



ref. [12 



In summary, our main finding is as follows: Improving the vertex Ansatz does little to 
improve the analytic structure of the quark propagator. However, improving the Ansatz 
employed for the gluon propagator, particularly by including a realistic asymptotic ultravi- 
olet tail, moves poles in the heavy quark propagator to a less intrusive part of the complex 
momentum plane. This portends well for future application of the heavy quark DSE tech- 
nique. 

The layout of the paper is as follows. In Section II we summarise the renormalised 
quark DSE and the approximations we shall be employing for the quark-gluon vertex and 
gluon propagator. In Section III we summarise the heavy quark formalism and derive a 
renormalised DSE for the heavy quark propagator. Numerical solutions to these equations 
are discussed in Section IV. Conclusions are drawn and suggestions for the direction of future 
work are given in Section V. 

II. THE QUARK DYSON-SCHWINGER EQUATION 

Our starting point is the renormalised quark DSE |1|] 

An 2 rA fj 4 a 
E'(p,A) = Z 1 ( / i 2 ,A 2 )^-y —S-D^p-q^S^T^p), (2.1) 

where we have used a Euclidean metric in which timelike vectors satisfy p 2 = — ^Minkowski < ^' 
and for which {7^,7^} = 2S^„. Our aim is to solve the DSE for the renormalised quark 
propagator S(p,fi), which we write in the form 

s(p, n) 



ry • pA(p 2 , /i 2 ) + B(p 2 , /1 2 ) 
1 



Z 2 ( / x2 )A 2)[, 7 . p + mo (A)] + S '( P)A )- 
The unrenormalised self energy is written 



(2.2) 



E'(p,A) = i 7 .p A'(p 2 ,A 2 )-l +B'(p 2 ,A 2 ). (2.3) 



If the renormalisation scale is set such that 



S(p)\ p i = s = , -t 1 , 2V (2.4) 

p M 27 • p + m R (n z ) 

it follows from Eqs. (|2.2j) and (|2.3|) that renormalised and bare quantities are related by 



Z 2 (/AA 2 ) = 2-AV,A 2 



(2.5) 



a(pV) = i + aV,a 2 )-^V,a 2 ), 



(2.6) 



B(p 2 ^ 2 ) = m R (^) + B\p\ A 2 ) - B'(^ 2 ,A 2 ). 



(2.7) 



The set of equations ( |2.1| ) and ( |2.5| ) to fl2.7| ) together with the 'abelian approximation' 
Z\ = Z2 can be solved numerically for the propagator functions A and B once the renor- 
malised quark-gluon vertex function T^, the renormalised gluon propagator D^, the renor- 
malisation point (/i, J71r(^)) and cutoff A are specified. The precise forms of the quark-gluon 
vertex and gluon propagator are unknown and must be modelled by appropriate Ansatze. 
We next summarise the Ansatze employed in this paper. 

A. Quark-gluon vertex 

The most general form of the quark-gluon vertex consistent with Lorentz and CPT 
invariance, satisfying the Ward identity iY^p^p) = d 1 ^S~ 1 {p) and Ward Takahashi identity^] 
i(p — g) M r M (p, q) = S^ 1 (p) — S~ l (q), and free of kinematic singularities has been given by 
Ball and Chiu 111 511 . It takes the form 



r M (p,g)=r? c (p,g)+rJ(p,g), 



(2.8) 



where 



1 r 



rrM = i Mp 2 ) + A( q 



+ 



(p + g)/ 

(p 2 — q 2 



7^ 



A{p 2 ) - A{q 2 



1-P + l-Q 



B{p 2 ) - B{q 2 



(2.9) 



and rT(p, q) is an otherwise unconstrained piece satisfying the conditions (p— g) M rJ(p, q) = 



and rj(p,p) = 0. 

We mention two well studied vertex Ansatze falling within this class. The first of these, 



introduced by Curtis and Pennington jl7| to ensure multiplicative renormalisability in quan- 



tum electrodynamics, is defined by setting the transverse piece T^ equal to 



rj cp (p,g) 



A( P 2 ) - A(q 2 ) 
2d(p,q) 



7/i(p 2 - <f) - (p + q)p(l -p-l-q) 



(2.10) 



with 



1 Strictly speaking, it is the Slavnov- Taylor identities, which include ghost contributions, and not 
the Ward- Takahashi identities which are relevant to QCD. By using the Ball-Chiu vertex we are 
effectively ignoring the ghost self energy and ghost-quark scattering kernel. This is a commonly 
used approximation M. 



d{p, q) 



{p 2 -q 2 ) 2 + [M 2 (p 2 ) + M 2 {q 
p 2 + q 2 



2/^2\21 2 



where M = B/A. 

The second of these, proposed by Haeri ||18| , takes the form 



p 2 — q 



p 2 — q 2 



B(p 2 )-B(q 2 
p 2 — q 2 



(7 • VI i, + 7 M 7 • q) ■ 



(2.11) 



r^)- p2Aip v q : A(q2 \ ^ A{p v A y\ - 



pi^i • q 



(2.12) 



It satisfies the above criteria and therefore must be of the form Eq. (|2.8| ). Munczek flll| has 
shown that the Haeri vertex is identical to the spectral representation of the vertex used in 



the gauge technique, and this has in turn been employed in ref. [13 



B. Gluon propagator 



In a general covariant gauge, the gluon propagator takes the form 
g 2 D ilv {k) = I ,\, 



'{IV 



kfik ^ A{k 2 ) + g 2 ^ kv 



k 2 



k 4 



(2.13) 



where £ is the gauge fixing parameter. Perhaps the simplest Ansatz which has proved useful 
for modelling QCD is the 'infrared dominant' model 



(2.14) 



Am(fc 2 ) = -(27r)VV(fc). 
Id 



In applications to hadronic physics, \x is usually taken to be of the order of 1 GeV |20 



which is the typical scale of QCD. An obvious disadvantage of this model is that it ne- 
glects completely any intermediate or ultraviolet behaviour. One can go some way towards 
incorporating some intermediate energy structure by using the computationally convenient 
gaussian model 



A G (k 2 ) 



— (2tt) -rt 
16 a z ir z 



-k 2 /a 



(2.15) 



We note that the infrared dominant gluon model Air, together with the minimal Ball-Chiu 
vertex Y^ijp^q) = r™(p, g) defines precisely the light quark model considered in ref. [7], 
whereas the gaussian model Ag has been used in our earlier rainbow approximation studies 
of heavy quarks |I2 1 . Below we shall explore the effect on the analytic structure of the heavy 
quark propagator of combining either Air or Ag with the Ball-Chiu vertex. 

A more realistic model gluon propagator which goes some way toward modelling the 
asymptotically free ultraviolet behaviour of QCD (neglecting logarithmic corrections) has 
been proposed by Frank and Roberts [|16|]. It takes the form 



A FR (A; 2 



4vr 2 d 



A<K 2 m 2 t 5\k) + 



3 fe 2 /(4m2) 

T 2 



(2.16) 



where d = 12/(33 — 2Nf), Nf = 3 is the number of light quark flavours, and m t ~ 0.69 GeV is 
a parameter fitted to a range of calculated pion observables. In our numerical calculations, 
we find it more convenient to consider a gaussian smeared version of the FR propagator 
given by 



m 2 rl „ 1 — p fe2 /( 4m t) 

a^ir z k 1 



, GFR (k 2 ) = W^f-e-* 2 /* + ^ 2 d^—— . (2.17) 



This will enable comparisons to be made with the pure gaussian model A^. 

In the next section we consider the heavy quark limit of the quark Dyson-Schwinger 
equation. Our earlier analysis of this limit [[T2] employed the heavily damped gluon propa- 
gator Aq, and so was free from ultraviolet divergences. In order to deal with a more realistic 
model, such as Eq. ( |2.17| ), it is necessary to develop a properly renormalised version of the 
formalism. 



III. THE HEAVY QUARK LIMIT 

A. Heavy quark propagator 

In the absolute limit of heavy renormalised quark masses, m R (fi 2 ) —>■ oo, the dressed 
quark propagator Eq. Q2.2| ) is dominated by the bare form S^ ie — ij ■ p + m R . However it 
is important to isolate from the full inverse propagator finite order self energy corrections 
to the bare inverse propagator which drive confining and remnant chiral symmetry breaking 
effects. To this end we write the momentum variables occurring in the DSE Eq. (|2.1| ) as 

Pfi = im R (fi 2 )v^ + kfj,, q^ = im R (fi 2 )v^ + k'^ (3.18) 

where, for convenience, we take v = (0, 1), so k ■ v — k^. We then write the renormalised 
quark propagator functions as 

A(P 2 ,^) = 1 + ^A, (3.19) 

B(p 2 ,fi 2 )=m R (fi 2 ) + J: B (K,K), (3.20) 

where we have defined the independent momentum variable 

K = p2 + m « =h + J-, (3.21) 

2im R 2im R 



and renormalisation point 



V 2 + m R 
2im R 



(3.22) 



in the complex K plane. The change of dependent variable p 2 — > K induced by the transfor- 
mation Eq. ( |3.21| ) is illustrated in Fig. 0. In general, when working to zeroth order in l/m R , 



one can use the approximation K = k± (an approximation which was used in ref. [[L2|). An 
exception to this rule, relevant to models such as Eq. ( p. 17 ) for which the gluon propagator 



is not heavily ultraviolet damped, is in the denominator of the heavy quark propagator in 
the integrand in the quark DSE. This point will become clearer at the end of Section IIIB. 
From Eqs. (p^|), ( |3T9|) and (^) we have 

S(P, /*) = l V"- K Jr (K \ + O (—) , (3.23) 

2 iK + h(K,K) \m R J 

where we have defined the heavy quark self energy 

e(k,k) = z b (k,k)-e a (k,k). (3.24) 

We find in general that the DSE leads to a single integral equation for the complex valued 
function Yi(K, k). The form given by Eq. (|3.23| ) represents the heavy quark propagator in the 
dominant region near the the bare propagator mass pole p 2 = —m 2 ^. Obtaining an integral 
equation for 'E(K) involves the change of integration / d 4 q — > J d A k' = /f^ dk' 4 / °° d |k| |k| 
induced by the change of variable Eq. Q3.18| ). For this change of integration to be valid, the 
propagator, and hence the functions A(q 2 )/[q 2 A(q 2 ) 2 +B(q 2 ) 2 ] and B(q 2 )/[q 2 A(q 2 ) 2 +B(q 2 ) 2 }, 
must be analytic over the shaded region in Fig. pD, 

Re(q 2 )>-m 2 R + ^^. (3.25) 

Equivalently, the function defined by 

^ K '^ = wrwyy (3 ' 26) 

must be analytic over the shaded region ImK < 0. 

The confinement criterion that S(p) should be free from timelike poles on the negative 
real p 2 axis translates in the heavy quark case to a requirement that <tq should be free 
from poles on the imaginary K axis. The stronger conjecture 0, that the quark propagator 
should be an entire function of p 2 translates in the heavy quark formalism to a conjecture 
that <tq should be an entire function of K. 

From Eq. (|3.23|) , we have that the renormalisation condition Eq. (|2.4j) is equivalent to 



S(p,n)\ 



1 + 74 1 



"■ = " 2 2 iK 



(3.27) 

K=K 



to zeroth order in I/iur. Typically we choose k to be on the negative imaginary K axis, 
as the heavy quark propagator asymptotes to the bare propagator asiC^ — ioo, as can be 
seen from Fig. [I]. 

B. Heavy quark DSE 

To illustrate the derivation of the renormalised heavy quark DSE, we choose Landau 
gauge (£ = in Eq. (|2.13| )) and, for the time being, work with the rainbow or bare vertex 
approximation 



IV(p, q) = 7m- 



(3.28) 



Using Dirac trace identities to project out from Eqs. ( |2.1| ) to ( |2.3| ) a pair of coupled integral 
equations gives 



1 + 



A'(p 2 ,A 2 ) = 
AZ 1 r A d 4 q 



3p 2 J (2vr) 



p- q + 2 



p-{p-q)q-{p- q) 

(p-q) 2 



A 



(p-q)' 



A(q 2 



q 2 A 2 + B 2 '' 



(3.29) 



Substituting 



B'(p 2 ,A 2 )=AZ l [ A £^A 



(2ir)< 



(p-q)' 



B(q 2 



A(q 2 



q 2 A(q 2 ) 2 + B(q 2 ) 2 2m R iK' + Y,(K') 



q 2 A 2 + B 2 ' 



0\Ar 



m 



R , 



^3.30) 



(3.31) 



B(q 2 



1 



1 



q 2 A(q 2 ) 2 + B(q 2 ) 2 
into Eqs. (gjg) and ( ^30|) gives 

2Zi /" A d 4 A;' 



2 iK' + E(AT') + UijJ' 



m R 



A'(p 2 ,A 2 ) 



xA 



(2tt)< 



/\2 



1 + 2 



(k - k') 



(h - K) 2 

(k - k') 2 
1 



iK' + E(K', k) + \m R 



and 



B '^ A2 )- 2Z '/ A |^ A 



A2 



(k - k') 



iK' + E(K',K) + \m R ) 



From Eqs. Q, Q, (gl9| ), ( ^20|) and flOg ) we have 

S(ir, K )=[ J B'(p 2 ) A 2 )-m R A'(p 2 ,A 2 
which gives 



I? 2 - a 2 ] 



v ' ; 3 7 (2tt) 4 



|k-k'| 
(A; - k') 



:A 



/\2 



(A; - A/) 



iK' + ^(K',n) 



(K->k) 



(3.32) 



(3.33) 



(3.34) 



(3.35) 



(3.36) 



where i^ on the left hand side is related to k under the integrand via Eq. ( [3.2 1| ), and a 
similar definition exists relating K' and k' . Here we have assumed that Eqs. ( |2.5| ) and ( |3.33| ) 
imply Zi = 1 + 0(1 /m R ). 

As noted above, K and £4 are interchangeable to leading order, except where they occur 
in the denominator of the heavy fermion propagator l/[iK' + Jl(K')] in the integrand of 
Eq. (|3.36|) . This is because, for any gluon propagator A with a realistic asymptotic UV 



8 



behaviour, all powers of k' must be retained in the denominator to maintain the same 
degree of divergence in the heavy fermion DSE as in the original equation ( j2.1|) . With this 
observation, and choosing k^ = (0, K) and the renormalisation point (0, k), we arrive at the 
leading order heavy quark DSE 

4 M d A k' 1 



n.K, K) = -f 



x 



(2?r) 4 iki + k' 2 /(2m R ) + E(*4, k 
{K - k\f + |k' 



IkfA 



(K-k' A y + \k 



/I 2 



(K^k)}. (3.37) 



We show in the Appendix that, with the smeared FR gluon propagator Ansatz Eq. ( |2.17| ), 
and assuming a hierarchy of scales 

m t , \K\ , |«| < m R < A, (3.38) 

the integral in Eq. ( |3.37| ) is independent of the ultraviolet cutoff A, and the heavy quark self 
energy behaves like 

E(K, k) ~ 2id(K - K) In (—Y (3.39) 

as m R — >■ oo. With the more severely truncated gaussian gluon propagator Eq. fl2.15|) , the 
renormalisation point K can be taken to — ioo and the k' /(2m R ) term in the denominator 
of the integrand ignored with impunity. 



C. Choice of renormalisation point 

At the end of the day, physical quantities must be insensitive to the choice of renormali- 
sation point K. In this section we note that the freedom to choose the renormalisation point 
is equivalent to the notion of a 'residual mass' in HQET |2l| , |i4"|| , that is, the notion that to 
zeroth order in l/m R , physical quantities computed in HQET do not depend on the choice 
of m R . 

After formally carrying out the spatial momentum integration in Eq. ( |3.37|) , one obtains 
an equation generically of the form 



1 r 00 

E(K, k) = — -iK= dK' [T(K - K') - T{k - K')} a Q (K', k), (3.40) 



for some kernel T, and with oq defined by Eq (|3.26| ). It is possible to show from this generic 



form that the effect of making a change of renormalisation point /t id — > Knew is equivalent 
to a shift of the quark propagator solution along the imaginary K axis: 

cr Q (K, K new ) = a Q (K - i5m, K \ d ), (3.41) 

where 5m is the solution to 

= 0-q(kW - iSm, K old ), (3.42) 

9 



that is, the shift is that required to ensure that the new heavy quark propagator pass 
through the renormalisation point cr<2(K new , K new ) = l/m new . Referring to Fig. we see that 



this is equivalent to a shift in the position of the origin of the K plane along the p axis 
corresponding to changing tur by an amount 8m. 

One can also demonstrate that the mass differences between any two heavy quark-light 
antiquark meson states calculated from the Bethe-Salpeter formalism set out in ref. [12| is 
independent of the renormalisation point. 



IV. RESULTS 



Our main concern in this paper is to compare how the analytic structure of the heavy 
quark propagator solutions is affected by the approximations employed both for the quark- 
gluon vertex and the gluon propagator. In an earlier work jl2| an attempt was made to study 
the heavy quark-light antiquark meson spectrum using a combination of rainbow DSE and 
ladder BSE. It was found that, if a simple gaussian Ansatz is used for the gluon propagator, 
complex conjugate poles occur in the heavy quark propagator which prevent solution of the 
meson BSE. Below we systematically explore the movement of the poles as the bare vertex 
of the rainbow approximation is replaced by a Ball-Chiu vertex, and as the gaussian gluon 
propagator is replaced by the more realistic Frank and Roberts propagator. 



A. Gaussian gluon propagator 



1. Rainbow approximation 



As noted in the previous section, if the gluon propagator Eq. (|2.15|) is used, we may set the 
renormalisation point k = —ioo and ignore the k' /(2m#) term in Eq. ( |3.37| ). Furthermore, 
there is no need to distinguish between the independent momentum variables K and k^. 
With these simplifications the Landau gauge, rainbow heavy quark DSE becomes 



S(fc 4 



4 r ofk! |k-k'| 2 A((k-k') 2 ) 
3 J (2tt) 4 (k - k') 2 ik' A + £(A4) ' 



(4.43) 



Choosing the infrared dominant gluon propagator Eq. fl2.14| ), the DSE reduces to an algebraic 
equation with solution 



S(fc 4 



— ik^ + 



*-f-kl) if 0<A;4< , 



Vsp 



-*vHMi-¥) if fc 4 ># 



2 ' 



or, using the definition ( |3.26|) 



oq(>4 



3^ 2 



8i 
3fi' 2 



—iki + 



:V 



- kl) if < k A < 



v^M 



-h + \lkl - % 



if k 4 > 



2 ' 



(4.44) 



(4.45) 



10 



Alternatively, choosing the gaussian gluon propagator Eq. (|2.15| ) and carrying out the 
:grat: 

£(fc 4 ) 



c? 3 k integration, we obtain the integral equation [22 

,2 



/' 



2a 2 y/n 

-(h- 



dk A 



J\2 



2/v 4 T 2j^fc 4 

-(k 4 -k' 4 f/a 



X 



*> /a + ^\k A -k' A \ei{c 



i |3 __r_ / 1^4 — fe 4 | 



O 



(4.46) 



where erfc z = 1 — erf z is the complementary error function. This equation can be solved 
numerically. 

2. Ball- Chiu vertex 



If any of the minimal Ball-Chiu Ansatz Eq. (|2.9|) , the Curtis- Pennington Ansatz Eqs. 
and ( |2.10| ), or the Haeri Ansatz Eq. ( ^.12[ ) is used in place of the bare vertex, together with 
the Landau gauge gluon propagator, we obtain in place of Eq. ( |4.43|) the equation 



if ^L |k - kf A ((£; - A;') 2 ) 
( A > ~ 3 J (2vr) 4 (A; - k 1 ) 2 ik' A + S(fc 4 ) 



1 + 



S(fc 4 )-S(A; 4 ; 
i(k A — k' A ) 



(4.47) 



It is interesting to note that, within the set of vertex Ansatze we have considered, the heavy 
quark propagator is insensitive to the transverse part of the vertex. This is not difficult to 
understand for the Curtis-Pennington vertex, in which the transverse part is heavily damped 
by the presence of the factor M 4 ~ m 4 in the denominator d(p, q). However, in the case of 
the Haeri vertex there is no such obvious mechanism, and one is led to question whether 
the heavy quark propagator may be insensitive to a broad class of Ansatze satisfying the 
criteria specified above Eq. Q2.8p. 

Taking the gluon propagator to be the infrared dominant form Eq. (|2.14j) , gives the 
differential equation 



E(fc 4 



3/i 2 d 
16z dki 



\n[ik A + S(/c 4 )], 



(4.48) 



which, together with the boundary condition o"q(/c 4 ) —* as A; 4 —* —ioo, admits the solution 



a Q (k A ) = (3 V^e' 13 fc * + 2iF(-pk 



where gq is defined by Eq. ( p.26|) , (3 = 2v / 2//Uv / 3 and 



F(z) = e~ z2 [ e t2 dt = l -^L e - z2 erf{-iz) 
Jo 2 



(4.49) 



(4.50) 



is Dawson's integral. We note that this solution is an entire function of fc 4 , which, as pointed 
out earlier, is a desirable feature of a quark propagator. This comes as no surprise, as it 
simply the heavy quark limit of the model considered in ref. J7J, in which it was demonstrated 
that the combination of Ball-Chiu vertex and infrared dominant gluon propagator leads to 
an entire function propagator for all values of the bare current quark mass. 
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It is of interest to determine to what extent this analytic structure is a feature of the Ball- 
Chiu vertex, and to what extent it is a feature of the infrared dominant gluon propagator. 
If the infrared dominant propagator is replaced by the gaussian smeared form Eq. ( [2.15 ), 
we obtain the integral equation 



E(fc 4 



F 



2a 2 0F 



dk/i 



1 



2/C4 ~t 2-i\K^ 



X 



(k 4 



Kf 



-{k^k' A f/ a 



+ \fn\k A 



k' 4 1 erfc 




x 



1 + 



s(fc 4 ) - E(fc 4 ; 

i(h - k' 4 ) 



(4.51) 



which can be solved numerically. 



3. Numerical results: gaussian gluon propagator 



For the purpose of determining the analytic structure of the heavy fermion propagator 
obtained from the DSE with the gaussian gluon propagator Eq. ( |2.15| ), it is sufficient to look 
at the one parameter family of models obtained by scaling either /1 or a to unity. We choose 
to scale jjl to unity, which amounts to working with a set of dimensionless quantities 



&4 = k^/fi, 



a 



a /n 2 , &Q = 1*7q. 



(4.52) 



This choice enables us to recover the infrared dominant model in the limit a — > 0. 

In Figs. and |^ we plot the heavy quark self energy S(/c 4 ) as a function of real £; 4 
obtained from the bare vertex DSE Eq. Q4.46|) and the Ball-Chiu vertex DSE Eq. (|4.51|) for 
a — 1, 2 and 3. These results are obtained by iterating from an initial guess and using 
a Simpson's rule quadrature. We find that the derivative-like terms in Eq. ( |4.51|) prevent 
a numerical solution for values of a less than 1, as numerical noise in the function values 
becomes unstable with respect to iteration at small values of fc 4 . This problem is a general 
feature of numerical treatments of DSEs with Ball-Chiu-like vertices. Also plotted are the 
a = analytic results Eqs. ( [4.44 ) and ( (4.49| ). In all cases the self energy is characterised by 
a real part which peaks at zero and an imaginary part which peaks near the typical scale of 
the model /c 4 ~ fi. The self energy for negative real fc 4 can be obtained from these results 
using the reflection property £(— fc|) = S(/c 4 )*. 

To solve for the heavy quark propagator away from the real fc 4 axis we shift the contour of 
integration into the complex plane parallel to the real fc 4 axis and again solve iteratively. We 
note that, to determine S(/c 4 ) for complex arguments, it is necessary to move the contour 
of integration to pass through the point /c 4 . This is because the radial part of the d 3 k.' 
integration, carried out in going from Eq. (|4.43|) to Eq. (|4.46|) or from Eq. ( (4.47|) to Eq. Q4.51D , 
creates a pinch singularity at k' 4 = k 4 in the error function term in Eq. ( |4.4(j| ) or Eq. ( 4.51 ). 

We have carried out a search for poles in the propagator function <tq(/c 4 ) for a range 
of values of a for both the bare and Ball-Chiu vertex. Our results are listed in Table |. 
In all cases we find that the only observed poles occur for Im/c 4 > 0, and that gq dies 
away to small values and is free from singularities over that part of the shaded region in 
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Fig. m accessible to our computer program. Of course we are unable to pass the contour of 
numerical integration through the pole itself, and these results are attained by extrapolation 
from results of contours which we gradually moved deeper into the complex plane. We were 
unable to obtain a reasonable extrapolation for the Ball-Chiu vertex at a = 1, again because 
of the iterative instability problem associated with the derivative-like term in Eq. ( [4.51 ) . 
We also list in Table | the results of using the model gluon propagator 

g 2 D^(k 2 )=5^A G (k 2 ), (4.53) 



where Aq is given by Eq. (2.15). Forms such as Eq. ( |4.53p are frequently used in phenomeno- 



logical modelling (see for instance ref. 0) and are sometimes referred to as propagators in 
a 'Feynman-like gauge', though of course they are generally not of the form of Eq. Q2.13|) . 
The sole advantage of the Feynman-like gauge is that it leads to considerably simplified 
calculations. In our case it is possible to locate poles more accurately because there is no 
pinch singularity requiring the contour of integration to pass through the point in question. 
Once the propagator has been solved on the real k^ axis, the value of the propagator can be 
calculated at any point in the complex plane by integrating once along the real axis. Never- 
theless, we have also repeated our pole calculations by shifting the contour and extrapolating 
as in the Landau gauge case as a check on the the consistency of the two methods and find 
that they agree to within the accuracy given in Table | of the corresponding Landau gauge 
results. A Feynman-like gauge propagator was also used in ref. [|T^] dealing with the ladder 
Bethe-Salpeter equation for the heavy quark-light antiquark system. There it was demon- 
strated that the model with bare quark gluon vertex and gaussian Feynman-like gauge gluon 
propagator had no solutions because of poles in the heavy and light quark propagators. Ide- 
ally, one would like improvements in the DSE approximations to move the poles further from 
the real k^ axis to avoid the region of the k^ plane sampled by a Bethe-Salpeter calculation. 
From Table | we conclude that simply replacing the bare vertex by the Ball-Chiu vertex 
in itself does nothing to improve the pole structure of the heavy quark propagator, either for 
the Landau gauge gluon propagator or the Feynman-like gluon propagator. In particular we 
find that, as the gaussian width a increases, a mass pole pole moves in along the imaginary 
&4 axis. A pole on the imaginary axis indicates that the fermion can propagate as a free 
particle, and the position on the positive imaginary axis gives the contribution to the quark 
mass from the dynamical self dressing. For the bare vertex, the pole splits into conjugate 
pairs either side of the imaginary axis as a decreases. In this instance the quark becomes a 
confined particle. Numerical difficulties described above prevented us from confirming that 
the same situation occurs in the case of the Ball-Chiu vertex. As a — > we must recover 
the solution Eq. ( 4.49| ), which is an entire function with an essential singularity at infinity. 



We see from Table [I] that, when poles occur, their position remains almost unchanged 
in going from Landau to Feynman-like gauge if the Ball-Chiu vertex is used, but not if the 
bare vertex is used. In a properly formulated gauge covariant calculation, the position of 



any propagator mass pole should be independent of the gauge fixing procedure p3fl . While 
we certainly do not claim that that our treatment is gauge covariant, it is amusing to 
note that replacing the bare vertex by the Ball-Chiu vertex (and Feynman gauge by the 
computationally convenient Feynman-like gauge) appears to go some way towards satisfying 
this requirement. 
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B. Prank and Roberts gluon propagator 

We now return to the renormalised DSE Eq. (|3.37j ) with the gluon propagator A set equal 
to the smeared FR Ansatz Eq. ( [2.171) . For numerical simplicity we shall restrict ourselves 
to the rainbow approximation r M (p, q) = 7 M - We set 

K = X + iY, K = ir). (4.54) 

Assuming the contour of integration can be deformed to pass through k' 4 = K for the integral 
of the first term in chain brackets in Eq.( |3.37| ), and through k' A = k for the second term, we 
further set 

k' 4 = x + iY and k' 4 = x + ir], (4.55) 

respectively in each of these two terms. We also make the replacement k' 2 /(2m R ) — > (x 2 + 
|k'| )/(2mji) without affecting £ to leading order. This gives 

dx r d 3 k' 



E(X + iY,i V ) = U°°^-[ 

6 J-oo Zti J 



(2tt) 3 



1 Ik'l 2 



ix - Y + {x 2 + \k'\ 2 )/(2m R ) + E(x + iY) (x - X) 2 + |k'| 2 



A 



[x-X) 2 + \k'\ 



Ikf 



ix — rj + (x 2 + |k'| )/(2m R ) + E(x + irf) x 2 + |k'| 



A 



x 2 + |kf . (4.56) 



For the purpose of carrying out the numerics, it is convenient to change to the polar 
coordinates x = rcos0, |k'| = rsin0, giving finally 



I /"OO f"K 

T,(X + iY,ir)) = — - dr d(f)r 3 sin 
3tt 6 Jo Jo 



*■■- 2 (f> 
1 r 2 sin 2 . 



ir cos (f) — Y + r 2 /(2rrifj) + £(r cos cj) + iY) r 2 — 2Xr cos cj) + X 2 
1 



A(r 2 -2Xrcos0 + A 2 ) 



ir cos — r] + r 2 /(2m R ) + E(r cos + ir]) 



sin 2 0A(r 2 ) } . (4.57) 



This equation is first solved numerically along the line \mK = rj, (i.e. Y — 77), and the 
function E along this line is stored for subsequent calculations at arbitrary Y. 

1. Numerical results: Frank and Roberts propagator 

We have numerically solved Eq. (|4.57|) over a region of the complex K plane with the 
smeared FR gluon propagator Agfr- Our parameter choices are m t = 0.69 GeV, in agree- 
ment with ref. [16|, and a = 16m 2 ad/3 = 0.5643 (GeV) 2 corresponding to a — 0.5. The 
choice of a is designed so that a comparison can be made between the full FR propagator 
and the gaussian propagator Eq. Q2.15 ) obtained by keeping only the first term in Eq. ( p. 17 ). 



From Table | we know that gaussian propagator Aq with the parameter choice a = 0.5 re- 
sults in a pair of conjugate poles in the heavy quark propagator. The heavy quark mass is 
set to m R = 5.0 GeV. 
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In Fig. f| we plot the modulus \ctq(K) | of the heavy quark propagator for the full gaussian 
FR propagator A GF r and in Fig. || plot the same quantity using only the gaussian part of 
Aq with the parameters otherwise unchanged. Both calculations have been done using the 
renormalisation point k = — l.Oi GeV, and to clarify the comparison the same region of the 
K plane is displayed in both plots. The calculation of <jq(K) involves a shift of integration 
path to a contour parallel to the real axis passing through the point K. Deforming the 
contour to include points behind the pole is a numerically tedious exercise which is unlikely 
to enhance our understanding, so no results are given for the part of the K plane in Fig. |5|. 

In the process of carrying out our computations, we have observed that the shift property 
resulting from changes of renormalisation point, namely Eq. ( |3.41| ), is indeed respected by 
our numerical solutions. In fact, the full plot in Fig. f| was pieced together by altering the 
renormalisation point to obtain solutions in strips of the complex plane parallel to the real 
axis, and using Eq. (|3.42| ) to match solutions where strips overlapped. 



We note a clear movement of the propagator pole further away from the real K axis when 
the asymptotic ultraviolet tail is included in the gluon propagator. Since it is the proximity 
of the pole to the real axis which prevented a solution to BSE in our earlier studies, this 
suggests that this movement of the pole portends well for future possible studies of heavy 
quark mesons if careful attention is paid to the asymptotic ultraviolet behaviour of the gluon 
propagator. 

However, a note of caution is in order. A different choice of renormalisation point would 
result in the plots in Figs. [| and || shifting by different amounts respectively along the 
imaginary axis, ensuring that both plots pass through the same point ctq(k) = 1/ik. Conse- 
quently, the actual amount by which the pole moves away from the real K axis as a result of 
adding an asymptotic tail to the gluon propagator is an artefact of the choice of renormali- 
sation point, though the movement will always be away from the real axis. Of course only 
by carrying though the BSE calculation completely can one say for certain whether bound 
state meson solutions can be obtained. 

V. CONCLUSIONS AND OUTLOOK 

We have explored the analytic structure of heavy quark propagators following a recently 
proposed formalism which borrows ideas both from the DSE technique and HQET. It is our 
belief that, if the successes of HQET are to be properly understood, we must first understand 
how the non-perturbative dynamics of QCD affect the heavy quark propagator. Within the 
light quark sector the analytic structure of the quark propagator is perhaps best understood 
in terms of model Dyson- Schwinger equations. It is therefore a worthwhile exercise to extend 
the DSE technique to the heavy quark limit. 

The initial attempt in this direction [Q failed essentially because the approximations 
used led to spurious propagator poles which prevented solution of the bound state Bethe- 
Salpeter equations. Two approximations were involved: modelling of the quark-gluon vertex 
and of the the gluon propagator. We have focused on each of these aspects in turn in this 
paper. In order to deal with an improved Ansatz for the gluon propagator with a realistic 
asymptotic ultraviolet behaviour, it has been necessary to formulate a properly renormalised 



version of the heavy quark DSE technique proposed in ref. [12]. As an interesting corollary to 
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our formalism we observe that the freedom to choose the renormalisation point is tantamount 

to the freedom in zeroth order HQET to choose the heavy quark mass up to a residual mass. 

We have first examined the effect of replacing the bare vertex with Ansatze based on 



the Ball-Chiu vertex | I5| , which is primarily designed to satisfy the Ward-Takahashi iden- 



tity. Specifically, we have considered the minimal Ball-Chiu vertex and two variants: that 



proposed by Curtis and Pennington |L7[ and that proposed by Haeri [|n|. The two variants 
differ from the minimal vertex by the inclusion of extra transverse components. We find 
that, to zeroth order in the inverse of the heavy quark mass, the heavy quark propagator 
is insensitive to which of the three above Ansatze is used. One is led to question to what 
extent transverse additions to the minimal Ball-Chiu vertex can be ignored in determining 
the leading order heavy quark propagator. If they can be be ignored in general there are 
immediate benefits in using the heavy fermion limit as a test-case for studies of confining 
field theories. 

In our numerical calculations we began with the Landau gauge form of the model gaus- 
sian gluon propagator A^, Eq. (|2.15|) , which was employed (together with the bare vertex) in 
previous studies ||12| . Unfortunately, we find no improvement in the propagator pole struc- 
ture in going from the bare vertex to the Ball-Chiu vertex while maintaining a gaussian 
gluon propagator. That is to say, timelike mass poles indicating non-confinement, or conju- 
gate poles which are likely to interfere with the successful solution to bound state problems, 
are not removed simply by improving the quark gluon vertex Ansatz alone. However, in 
the limit in which the width of the gaussian gluon propagator is taken to zero (the 'infrared 
dominant model'), we do obtain an entire function heavy quark propagator, free from sin- 
gularities except an essential singularity at infinity. This is consistent with the equivalent 
finite quark mass calculation |jj , and may provide a useful propagator for phenomenological 
modelling purposes. 

We conclude then that it is most likely the remaining approximation, namely the gaussian 
model gluon propagator, which is responsible for the poor analytic structure previously 
obtained for the quark propagator. To explore this possibility, we have replaced the simple 
gaussian gluon propagator Ansatz of ref. [Q by a gaussian smeared version of the more 
sophisticated Frank and Roberts Ansatz Apr given by Eq. fl2.16| ). In this case, convergence 
of the integral in the DSE which is lost by naively retaining only the lowest order of the 
1/rriR expansion of the quark propagator must be restored by judiciously including at least 
the spatially dependent 0(l/m/j) part: 

SOp) = i^i * + o (rrfZ 1 ) . (5.58) 

The renormalised current quark mass ttir then becomes an ultraviolet regulator, and in 
the case of the FR propagator the lowest order contribution to the mass expansion of the 
heavy quark self energy behaves as ln(mii/mt), where m t delineates the scale at which the 
asymptotic ultraviolet behaviour of the FR propagator sets in. 

Our numerical solutions of the heavy quark DSE show a clear movement of the offensive 
propagator poles away from that part of the complex momentum plane likely to be sampled 
by a Bethe-Salpeter calculation of heavy quark meson states. However we caution that the 
amount by which the poles shift is, strictly speaking, dependent on the choice of renormali- 
sation point. Without carrying through the the Bethe-Salpeter analysis one cannot say for 
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sure that the problem is solved. A further study of the BSE for heavy mesons is expected 
to be the focus of future work. 

Our work has also thrown up a couple of other interesting questions worthy of attention. 
Firstly, it should be possible to check directly to what extent the propagator pole structure 
is invariant with respect to the choice of gauge fixing parameter £. As the postions of 
propagator poles should be gauge independent p3| , this provides a straightforward measure 



of the ability of a particular vertex Ansatz to respect gauge covariance of the model. 

Secondly, one is led to question the meaning of propagator poles obtained from a Eu- 
clidean DSE formalism. Whether propagator poles obtained in this way are an artefact of 
the approximations used or whether they are a genuine property of quark propagators has 
been an open question for some time ||Tl,[T0|pl- It is possible that the heavy fermion limit may 
help to shed some light on this problem by developing a Bethe-Salpeter formalism for heavy 
quarkonium states. It is well known that the non-relativistic limit of the Bethe-Salpeter 
equation for a heavy fermion-heavy antifermion bound state can be written in the form of 



a Schrodinger equation [pq| . The derivation typically assumes physical mass poles in the 
fermion propagators whose residues contribute to the resulting Schrodinger equation. An 
analogous derivation for the case of propagator poles which have moved off the timelike 
momentum axis as a result of a confining gluon propagator may help both with interpreta- 
tion of quark propagator poles and with understanding the success of heavy quark potential 
models. 
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APPENDIX: CONVERGENCE OF THE HEAVY QUARK DSE 

Consider the heavy quark DSE Eq. ( |3.37| ), and suppose we assume for the gluon propa- 
gator Ansatz the asymptotic ultraviolet behaviour 

A(k 2 ) ~ ^jf for k 2 > ml (A.l) 

where m t is a scale parameter typically of the order of 1 GeV. This is the behaviour exhibited 
by the smeared FR Ansatz Eq. (|2.17| ). We demonstrate here that the right hand side of 
Eq. ( |3.37| ) is the finite difference of two logarithmically divergent integrals. 
We begin with the change of variables 

k' A = rcos<p, |k'|=rsin0, (A. 2) 



2 In an interesting recent development, McKay and Munczek p4l have examined the analytic 
structure of quark propagators when an extra constraint that solutions of the DSE should be 
Fourier transformable is imposed. 
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and hence 



d 4 k' = 4vr [dk' 4 fd\k'\ |kf = 4vr / dr 

J J Jo JO 



>r 3 sin 2 i 



(A.3) 



This gives 



Z(K,k) 



1 r A 



dr 



r 5 sin 4 ( 



37r 3 jo Jo ir cos 4> + r 2 / (2mji) + T,(r cos i 
~A(K 2 -2Kr cos (j) + r 2 ) 



K 2 — 2Kr cos 6 + r 2 



-(*->«) 



(A.4) 



Implicit in this equation is a hierarchy of scales given by Eq. ( |3.38|) . Simply counting powers 
of r in the integrand, we see that each of the two terms diverges as 



rAdr 



f " r 1 A 

m R / m R mA. 

J r 



(A.5) 



On the other hand consider the difference of the two terms. For r > m t the part in 
square brackets is 



4vr 2 rf 



(K 2 -2Kr cos (f) + r 2 ) 2 
= 4n 2 d 



(K->k) 



A(K-k)cos(J) 2(K 2 -K 2 )(6cos 2 (p-l) 



+ o 4 ■ 



(A.6) 



Neglecting E(rcos0) in Eq. ( A.4 ) for large r, we can approximate the propagator contribu- 
tion to the integrand by 



1 icoscj) — r/2mR 

ir cos <ft + r 2 / (2m ft) r(cos 2 <fi + r 2 /4m|.) 



(A.7) 



Then, taking into account the hierarchy (|3.38|) , the contribution to the integrand for r > m t 
is approximately 



16id(n — K) r 00 dr r 71 sin 4 cos 2 . 



3tt 



m t r Jo cos 2 (f) + r 2 /Am 



(A.8) 



giving 

T,(K, k) ~ 2id(n — K) In ( ) as rrifi — > oo. 

\m t ) 

The last step can be achieved using the crude approximation 



cos 



cos 2 + r 2 /4m 2 



1 if r < 4m ft 

(4m 2 R /r 2 ) cos 2 (f) if r > 4m^ 



(A.9) 



(A.10) 
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TABLES 

TABLE I. Position k/±j \x of poles in the heavy quark propagator closest to the real k& axis, using 
the gaussian gluon propagator Ansatz Aq. Numerical instabilities prevent an accurate location of 
the pole in the case indicated by a question mark. 





Landau 


gauge 


Feynman-like 


gauge 


a = a/ /j, 2 


bare vertex 


BC vertex 


bare vertex 


BC vertex 


0.5 


±0.510 + 0.415i 


? 


±0.5466 + 0.5109i 


? 


1 


±0.378 + 0.506i 


? 


±0.3383 + 0.6758i 


0.4429i 


2 


0.466i 


0.31i 


0.3844i 


0.3130i 


3 


0.307i 


0.252i 


0.2880i 


0.2554i 


4 


0.25L 


0.219z 


0.2408i 


0.2210i 
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FIGURES 

FIG. 1. The change of variables p 2 — ► K used to represent the heavy quark propagator in the 
region of the bare fermion mass pole p 2 = —m 2 ^. For the change of variables in the DSE to be 
valid, the quark propagator must be analytic over the shaded region. 

FIG. 2. The heavy quark self energy X(&4), from the DSE in Landau gauge with a bare 
quark-gluon vertex and gaussian gluon propagator with parameters: /x = 1 and a = (solid 
curve), 1 (long dashes), 2 (short dashes) and 3 (dotted curve). The upper curves are ReS and the 
lower curves ImS. 

FIG. 3. The same as Fig. ^, except with the Ball-Chiu quark-gluon vertex. 

FIG. 4. The modulus |o"q(ET)| of the heavy quark propagator in the complex K plane obtained 
by solving the heavy quark DSE in rainbow approximation with a gaussian smeared Frank and 
Roberts Ansatz Aqfr for the gluon propagator. Input parameter values are given in the text. 

FIG. 5. The same as Fig. |I|, except with the gluon propagator Ansatz replaced by only its 
Gaussian part Aq, and all parameter values otherwise unchanged. 
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